Strange and non-strange quark mass dependence of elastic light resonances from 
SU(3) Unitarized Chiral Perturbation Theory to one loop 
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We study the light quark mass dependence of the /o(600), k(800), p(770) and A"* (892) resonance 
parameters generated from elastic meson-meson scattering using unitarized one-loop Chiral Pertur- 
bation Theory. First, we show that it is possible to fit simultaneously all experimental scattering 
data up to 0.8-1 GeV together with lattice results on decay constants and scattering lengths up to a 
pion mass of 400 MeV, using chiral parameters compatible with existing determinations. Then, the 
strange and non-strange quark masses are varied from the chiral limit up to values of interest for 
lattice studies. In these amplitudes, the mass and width of the p(770) and K*{892) present a similar 
and smooth quark mass dependence. In contrast, both scalars present a similar non-analyticity at 
high quark masses. Nevertheless, the /o(600) dependence on the non-strange quark mass is stronger 
than for the /«(800) and the vectors. We also confirm the lattice assumption of quark mass indepen- 
dence of the vector two-meson coupling that, in contrast, is violated for scalars. As a consequence, 
vector widths are very well approximated by the KSRF relation, and their masses are shown to scale 
like their corresponding meson decay constants. 

PACS numbers: 14.40.-n 12.39.Fe 13. 75. Lb 



I. INTRODUCTION 

Although QCD is well established as the theory of 
strong interactions, the fact that its coupling becomes 
large at energies below 1-2 GeV keeps the hadronic realm 
beyond the reach of perturbative calculations. In that 
regime, lattice methods are a useful tool to calculate 
QCD observables, although the discretization involved 
in this technique introduces complications of its own, 
in particular related to chiral symmetry breaking and 
the implementation of realistic small masses for the light 
quarks. Despite the remarkable success of lattice stud- 
ies, results on light meson resonances are few and usually 
obtained at very large quark masses compared with their 
physical values P, [3| . This is particularly so for the light 
scalars, very relevant for nuclear attraction, but whose 
calculations are hindered bythe so called "disconnected 
diagrams". Very recently [3|, an alternative technique, 
based on Chiral Perturbation Theory (ChPT) and disper- 
sion relations, has been applied to calculate the depen- 
dence of the /o(600) (or "sigma") and p(770) resonances 
on the pion mass - in practice, the average mass of the 
u and d quarks. Now the starting parameters are phys- 
ical and resonances appear in amplitudes that describe 
real data on tttt scattering. The predicted dependence for 
the p{770) compares remarkably well with previous and 
later lattice predictions. For the scalar sigma it shows a 
non-analyticity that should be taken into account when 
extrapolating future lattice data to physical values. In 
this work we extend this study to include the strange 
quark mass within an SU(3) ChPT formalism. Our aim is 
threefold: first, to confirm previous results within a more 
general formalism. Second, to analyze the dependence on 
the average mass of the u and d quarks of the if * (892) 
and «;(800) strange resonances. The latter, despite being 
a scalar, and very similar to the /o(600), is much more 
feasible for lattice calculations within the next few years 



[J| due to its nonzero isospin and strangeness. Third, 
we also study the dependence of all the /o(600), k(800), 
p(770) and if* (892) parameters in terms of the strange 
quark mass. Finally, let us remark that the dependence 
of hadronic observables, meson masses in particular, is 
not only of relevance for lattice calculations, but also for 
anthropic considerations Q or the study of the cosmo- 
logical variability of fundamental constants @ . 

Thus, in the next two subsections we introduce very 
briefly the basic notation of ChPT, explain the relation 
between pseudoscalar meson and quark masses, and re- 
view the unitarization procedure. In section II, we show 
the fits to the existing experimental data on elastic scat- 
tering as well as to lattice results on pion and kaon 
masses, their decay constants, and scattering lengths on 
the highest isospin channels. Section III is devoted to 
the dependence of light resonance properties on the non- 
strange quark masses. In section IV we then study the 
dependence with the strange quark mass and in Section 
V we present our summary and conclusions. 



A. Chiral Perturbation Theory 

As is well known, pions, kaons and etas can be iden- 
tified with the Nambu-Goldstone bosons (NGB) of the 
spontaneous chiral symmetry breaking of QCD. If quarks 
were massless, so should be the NGB and they would be 
separated by a mass gap of the order of 1 GeV from other 
hadrons, thus becoming the only relevant QCD degrees 
of freedom at low energies. Of course, quarks are not 
massless, but the u, d and s flavors have a sufficiently 
light mass to be considered as a perturbation. It is thus 
possible to write a Low Energy Effective Lagrangian out 
of pion, kaon and eta fields, known as Chiral Perturba- 
tion Theory (ChPT)0]. This Lagrangian is built as the 
most general derivative and mass expansion that respects 
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the symmetries of QCD, particularly its chiral symmetry 
breaking pattern. Except for the leading order, fixed by 
symmetry and the scale of spontaneous symmetry break- 
ing, all terms in the Lagrangian are multiplied by a low 
energy constant (LEG) that contains the information on 
the underlying QCD dynamics and also renormalizes the 
loop diagrams with vertices from lower orders. In this 
way, pion, kaon and eta observables are obtained as a 
model independent expansion in powers of momenta and 
masses over the chiral scale 47r/o ~ 1.2 GeV, where /o is 
the pion decay constant in the chiral limit (as it is cus- 
tomary, for quantities at leading order in the quark mass 
expansion we will use the subscript). 

In particular, partial wave amplitudes for elastic 
meson-meson scattering are obtained within ChPT as an 
expansion 

t{s)^h{s)+U{s) + ...., hk = Oip'"'), (1) 

where p denotes either momenta or meson masses. Ac- 
tually, these partial waves carry definite isospin / and 
total angular momentum J, but we have momentarily 
suppressed these labels for clarity. As we have just com- 
mented, the leading order ^2(5) corresponds to the cur- 
rent algebra results and only depends on the scale of 
spontaneous chiral symmetry breaking /q. The next to 
leading order ^4(5) contains one- loop diagrams made of 



vertices from the lowest order Lagrangian, plus tree level 
diagrams of 0{p'^). Within the SU(3) formalism, these 
tree level diagrams are multiplied by LECs, denoted as 
Li, which are independent of masses or momenta, and 
have been determined from different experiments. In Ta- 
ble H] we provide several sets for the eight Li that appear 
in meson-meson scattering to one loop. Those with an 
"r" superscript carry a dependence on the regularization 
scale fj, [zl, customarily chosen at ^ = Mp. Of course, 
that scale dependence cancels in the calculation of phys- 
ical observables. The values in the second column come 
from the "Main Fit" of a Ku analysis to two loops Q, 
whereas those in the third column come from the same 
reference, but to one loop. Naively one would expect 
the LECs obtained in our unitarized one-loop fits to lie 
somewhere in between these two sets of values, since uni- 
tarization reproduces one of the most relevant numerical 
contributions from the two-loop calculation, namely the 
s-channel leading logs. As one of our main interests is 
ttK scattering and the K*{892) and k(800) resonances, 
we also provide the values obtained from a very rigorous 
treatment of Kit scattering lengths in terms the Roy- 
Steiner dispersion relations [9] . The rest of the columns 
correspond to unitarized ChPT fits that we will explain 
in Sect HI] 



TABLE I: 0{p*) chiral parameters (xlO^) evaluated at /i = Mp. The second and third 
columns come from the two and one loop analysis listed in 3, where L4 and Le were set 
equal to zero. The fourth column comes from a careful ttK dispersive analysis using the 
Roy-Steiner formalism. The lAMIII column is one of the sets obtained from an older fit 
with the coupled channel JAM (only statistical uncertainties are shown). The columns 
labeled Fit I and Fit II correspond to the simultaneous fit to experiment and lattice data 
performed in this work, which are described in Sec |lT] together with their uncertainties. 



LECs 


Ref.fS] Oip"-) Ref.[8| 0(p*) 


Ref.f9] 


lAM III 


Fit I 


Fit II 


LI 


0.53 ± 0.25 


0.46 


1.05 ± 0.12 


0.6 ± 0.09 


1.10 


0.74 


LI 


0.71 ± 0.27 


1.49 


1.32 ± 0.03 


1.22 ± 0.08 


1.11 


1.04 


Lz 


-2.72±1.12 


-3.18 


-4.53 ± 0.14 


-3.02±0.06 


-4.03 


-3.12 


LI 


(fixed) 


(fixed) 


0.53 ± 0.39 


(fixed) 


-0.06 


0.00 


LI 


0.91 ± 0.15 


1.46 


3.19 ± 2.40 


1.9 ± 0.03 


1.34 


1.26 


LI 


(fixed) 


(fixed) 




-0.07±0.20 


0.15 


-0.01 


Lv 


-0.32± 0.15 


-0.49 




-0.25±0.18 


-0.43 


-0.49 


LI 


0.62± 0.20 


1.00 




0.84 ± 0.23 


0.94 


1.06 


LI+2LI 


0.62± 0.20 


1.00 


3.66±1.52 


0.7 ± 0.46 


1.24 


1.04 


2LI - LI 


0.35±0.57 


-0.57 


0.78±0.24 


-0.02± 0.20 


1.09 


0.44 



In this work, we are interested in the quark mass de- 
pendence of the amplitudes, which appears in ChPT 
through Lagrangian terms that contain the quark mass 
matrix M — diag(TO, m, rris), that is treated as a per- 
turbation. Note that we work in the isospin limit rh ~ 
•niu = nrid = {m^^^ + mF^^^)/2. Chiral symmetry is ex- 
plicitly broken by these mass terms and the NGB acquire 



masses, which, at leading order, read 0]: 

Mo^ = 2mBo, 

= (m + m,)Bo, 

Ml, - ^(m + 2m,)i?o. (2) 
Let us recall that the constant Bq is defined from the 
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values in the chiral limit of the chiral condensate and the 
pion decay constant as fohows: Bf^ — — < Q\qq\Q >o //q ■ 
and thus it carries no quark mass dependence. To one 
loop, there are some corrections, and the physical meson 
masses now read: 



Ml 
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(4) 
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(5) 



327r2/,^ 



log- 



P — TT, K, rj. 



Note, however, that all the quark mass dependence 
always appears through the leading order masses M^p 
defined in Eq.(I2|). As a matter of fact, this also happens 
in the ChPT amplitudes, which means that studying the 
quark mass dependence, keeping Bq fixed, is nothing but 
studying the meson mass dependence. In practice, and in 
order to get rid of the Bq constant, we will recast all our 
results in terms of masses normalized to their physical 
values: 



^phys 



^"0 IT phys 



MIJ2 



m 



s phys 



M. 



K phys 



^O.phys/2' 



(6) 
(7) 



Note that, from now on, a quantity with a "phys" sub- 
script refers to the value of that quantity in the physical 
case. Thus, in this work we will change quark masses, 
that, using Eqs.Q, (|H) and ([S]), imply a change in meson 
masses, which are the ones appearing explicitly in the 
ChPT scattering amplitudes. 

There are many advantages in using meson masses as 
the variation parameter, since, contrary to quark masses 
that have a complicated and scale dependent definition 
on the QCD renormalization scheme, meson masses are 
observables, with no scale dependence and a straight- 
forward physical interpretation. Actually, many lattice 




s phys 



FIG. 1: Top: The ratio (M^/M^phya)/(m/mphys) Bottom: 
phys) • Within the range of variation 
of this work, a relative variation of a quark mass can be also 
understood as the same relative variation in the corresponding 
meson mass squared to within ~ 10% accuracy. The continu- 
ous and dashed lines correspond to Fit I and II sets of LECs 
given in Table U 



results are also recast in terms of pion or kaon mass vari- 
ations. Unfortunately the simple relations in Eqs.® and 
([7|) are exact only when written in terms of the leading 
order masses Mqp, not the observable ones. Neverthe- 
less, the one- loop corrections become numerically small 
when taking ratios so that, to a good degree of approx- 
imation, the reader still can think in terms of physical 
meson masses instead of their leading order values. Ac- 
tually, in FigH] we show that within the range of quark 
mass variations that we will consider in this work, the 
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naive, but intuitive, relations 



m 



Mi 



^phys 

rris 

phys 



phys 



Mi 



keeping ms = msphys, (8) 
keeping m = mphys (9) 



are a very good approximation ~ within less than 10% 
error- to the correct ratios in Eqs.dH) and ([7]), that we ac- 
tually use. To make our presentation of the results more 
intuitive we will give, when possible, our results both 
in terms of quark mass variation and the corresponding 
meson mass variation. 

At this point we have to address the question of how 
much we can vary the quark masses before our approach 
breaks down. First we want the pion always lighter 
than the kaon and eta since otherwise the elastic ap- 
proximation would make no sense for tttt or Ktt scatter- 
ing. Second, ChPT seems to work for masses as high 
as 500 MeV, since we already know that it provides a 
fairly good description of low energy Ktt scattering, even 
though Mk ~ 500 MeV. Thus, when changing the non 
strange quark mass, keeping fixed, we will show re- 
sults up to < 440 MeV but not beyond, since then 
Mk — 600 MeV. Equivalently, this means m/rhphys < 9. 
Concerning the strange quark variation with m fixed, we 
will consider 0.7 < ms/ws phys < 1-3, since M^ barely 
changes and 400 MeV < Mk < 585 MeV. This ensures 
that the rriK + rriri is not below the m^* mass so that 
we would need a coupled channel formalism. Of course, 
the closer to the estimated applicability limits the less 
reliable our formalism will be. 

The SU(3) TTTT and Ktt one-loop amplitudes were first 
calculated in [ll|, although for technical reasons ex- 
plained in [l^ needed for the implementation of exact 
unitarity later on, we use the expressions in the appendix 
of [l2|, but written in terms of all physical constants 
M.^, Mk, Mrj, fiT,fK,fr) as explained in [l^. For com- 
pleteness we show here the decay constant dependence 
on meson masses. 
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to the Li,L2 and L3 constants only contain derivatives 
and thus are somewhat less relevant for the quark mass 
dependence, but more relevant in terms of s dependence. 

Finally, let us remark that despite the fact that their 
effect is encoded in the LECs, the ChPT amplitudes, 
being an expansion, cannot describe resonances and their 
associated poles in the second Riemann sheet. Actually, 
resonances are usually identified with a saturation of the 
unitarity constraints, which for elastic partial waves of 
definite isospin / and angular momentum J read: 



(10) 



Of course, for tttt and Ktt elastic scattering the most 
relevant quark mass dependence comes via M^r, Mk 
and /tt, Jk (since etas only appear in loops). Conse- 
quently, the LECs that play the most important role 
are L4,L^,Le and Lg, since they appear in Lagrangian 
terms that contain explicitly powers of the quark mass 
matrix. In contrast, the Lagrangian terms proportional 



Imtijis)^a{s)\tij{s)f 



\tijis)\ < l/<j{s), (11) 



where cr(s) = 2k/y/s and k is the center of mass momen- 
tum. The above equations imply that the partial wave 
can be recast in terms of a single phase or "phase shift" : 



tuis) = exjpiS I J (s) sin Sij{s)/ a {s) 



(12) 



In this work we are only interested in the (/, J) — 
(0,0), (1,1) and (2,0) channels for tttt scattering and 
(/, J) = (1/2,0), (1/2,1) and (3/2,0) for ttK scattering. 
For simplicity we will drop the IJ subindex when dis- 
cussing general properties of elastic partial waves. 

Note, however, that the ChPT expansion Eq.(IT|), be- 
ing basically a polynomial in energy, violates the bound 
in Eq. pTj) as the energy increases and cannot generate 
poles. Still, ChPT satisfies elastic unitarity perturba- 
tively: 



Imt2(s) = 0, Inii4(s) = a|t2(s)|^ 



(13) 



But, of course, elastic unitarity can be badly violated if 
the ChPT series is extrapolated close to a resonance. For 
these reasons, the resonance region lies beyond the reach 
of standard ChPT. However, we will see next that ChPT 
can be used in an alternative way. 



B. Dispersion relations, unitarity and ChPT 

Instead of simply extrapolating its series to higher en- 
ergies, ChPT can be used to calculate the subtraction 
constants of a dispersion relation for the two-body am- 
plitude. These constants correspond to the values of the 
amplitude or its derivatives at a low energy point where 
the use of ChPT is well justified. The remaining infor- 
mation to build the amplitude comes from the strong 
constraints of analyticity and unitarity. 

First of all, it is straightforward to rewrite the strong 
non-linear elastic unitarity constraint given in Ea.(|lip. 
as follows 



Iml/t{s) = -cr(s). 



(14) 



This means that, from unitarity, we know exactly the 
imaginary part of 1/t in the elastic region. We are only 
left to determine the real part of 1/t. 

Concerning the analyticity constraints, for simplicity 
let us consider first the case of two identical particles. 
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as in TTTT scattering. Then, the analytic structure in the 
complex s plane is rather simple: it has a "right" or 
"physical" cut on the real axis from threshold to +00, 
and a "left cut" from —00 to s = 0. By means of the 
Cauchy Theorem, a dispersion relation provides the am- 
plitude anywhere inside the cut complex plane in terms 
of a weighted integral of its imaginary part over the cuts. 

In our case, instead of t we are interested in a disper- 
sion relation for 1/t since we know exactly its imaginary 
part in the elastic region thanks to Ea.([T^. For con- 
venience, and since t2 is real, instead of 1/t we define 
G = f^/t, that also has a right cut (RC) and a left cut 
{LG). Since scalar waves are known to have dynamical 
Adler zeros in the low energy region below threshold, we 
will also allow for a pole contribution PG in G{s). All in 
all, we can write a dispersion relation for G{s) as follows 



G{s) = G(0) +G"(0)s + iG"(0)s 



(15) 



s 

TT 



, Im G(s') 



LG{G) + PC. 



In the elastic region, unitarity in Eq.()ll|). together with 
Eg. dT^ . allow us to evaluate exactly ImG = —<yt\ — 
—Imti on the RC. Note the three 1/s' factors - called 
subtractions - that we have introduced to suppress the 
high energy part and in particular the inelastic contri- 
butions, so that the integrals are dominated by the low 
energy region. But once the integrals are dominated by 
the low energy, it is well justified to use ChPT inside the 
integrals and thus, for instance, the LC integral to one 
loop ChPT is given by LG{G) ~ LC{~ti) + .... 

The price to pay for the three subtractions is that ana- 
lyticity only determines the function up to a second order 
polynomial G{Q)+G' {Q)s+\G" {0)s'^ . However, note that 
its coefficients correspond to the values of the amplitude 
or its derivatives at s = 0, where ChPT can be safely 
applied. In particular, to one- loop, G(0) ~ ^2(0) — ^4(0), 
G'(0) ~ ^'2(0) - t'M and G"(0) = -^^'(0), since t'i{Q) 
vanishes. Let us neglect for the moment the pole con- 
tribution, which is of higher order and only numerically 
relevant below threshold. Then one finds that all contri- 
butions can be recast in terms of the leading ^2(5) and 
next to leading ^4(5) ChPT amplitudes. Finally, we ar- 
rive at the so-called Inverse Amplitude Method (lAM) 



t{s) 



t2{s) - ^4(5) 



(16) 



Remarkably, this simple equation ensures elastic unitar- 
ity, matches ChPT at low energies, and, using LECs com- 
patible with existing determinations, describes fairly well 
data up to somewhat less than 1 GeV, generating the 
/9, K*, a and k resonances as poles on the second Rie- 
mann sheet. It has been shown 15| that the scalars can 
actually be generated mimicking the LEC, tadpole and 
crossed channel diagrams by a cutoff of natural size, and 
thus it is said that scalars are "dynamically generated" 
from, essentially, meson- meson dynamics (meson loops). 



In contrast, to generate the vectors, a precise knowledge 
of the LECs is needed, namely, of the underlying, non 
meson-meson QCD dynamics. 

Here we will update this description of experimental 
data but furthermore we will simultaneously describe the 
existing lattice results for decay constants and some scat- 
tering lengths. 

The lAM equation above is just the one-loop result, 
but it can be easily and systematically extended to higher 
orders of ChPT or generalized within a coupled chan- 
nel formalism [l3, [H, EBl, generating also the ao(980), 
/o(980) and the octet cj). However note that there is no 
dispersive justification for the coupled channel approach 
formula |32| and that is the main reason, apart from sim- 
plicity, why we have restricted our analysis to the elastic 
case. 

For completeness, and even though it will be negligible 
except for very high masses near the applicability limits 
of our approach, let us now include the pole contribution 
PC ignored so far. Its contribution can be calculated 
explicitly from its residue (iTj and, to one-loop, we find 
a modified lAM (mlAM) formula: 



,mIAM 



t^ 



A 



ml AM 



t2 

ti{s2) 



<4 -I- A^IAM 

{s2-Sa){s-S2) [t'2{s2)-t'^{s2)] 



(17) 



-SA 



where sa is the position of the Adler zero in the s-plane, 
and S2 its LO approximation. The standard lAM is re- 
covered for yl^-f-^A/ _ which holds exactly for all par- 
tial waves except the scalar ones. Above, and in the 
usual lAM derivation A™^^*^ was neglected, since it 
formally yields a NNLO contribution and is numerically 
very small, except near the Adler zero, where it diverges. 
However, if ^'"■f^-'^ is neglected, the lAM Adler zero oc- 
curs at S2, correct only to LO, it is a double zero instead 
of a simple one, and a spurious pole of the amplitude 
appears close to the Adler zero. All of these caveats are 
removed with the mlAM, Eq. pT)) . The differences in 
the physical and resonance region between the lAM and 
the mlAM are less than 1%. However, as we will see, 
for large the a and k poles "split" into two virtual 
poles below threshold, one of them moving towards zero 
and approaching the Adler zero region, where the lAM 
fails. Thus, we will use for our calculations the ml AM, 
although it is only relevant for the mentioned second a 
and K poles, and only when they are very close to their 
corresponding Adler zeros. 

Finally, we want to comment on the unequal mass 
case, since we also want to describe Ktt elastic scat- 
tering. The main difference now is that the left cut 
extends from -co to s = (Afi — M2Y ^ and also that 
there is a circular cut, centered at s = with ra- 
dius \Ml — M||. Again their main contribution comes 
from a region where ChPT can be applied. This time, 
however, t2{s) has two zeros instead of one, S2± = 
i (^Ml + Ml ± 2y/4M^ - 7M2.M2 + 4M4) , and the 
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modification to the lAM reads: 



A 



ml AM 



{s - S2- 



S2^ 



- sa) 



SA 



(.If N I N , ^4(s2+)i£(S2+; 



r 



(18) 



Once again we note that this modification will be numer- 
ically negligible except in the close vicinity of the Adler 
zero. The poles of the resonances under study will only 
come close to that region for very high values of the quark 
masses, in the limit of applicability of ChPT and our ap- 
proach. 

Before describing our fits, we want to remark that, in 
the lAM derivation above, ChPT does not play any role 
outside its applicability limits. By including three sub- 
tractions we have suppressed strongly all contributions to 
the integrals in high energy regions where ChPT results 
are not reliable. Finally, the three subtraction constants, 
which correspond to values of the amplitudes or their 
derivatives at s = are well calculated with ChPT. Of 
course, this is just a one-loop calculation, although the 
generalization to higher orders is tedious but straightfor- 
ward. Hence, our approach does not model the left or 
inelastic cuts, but just uses the corresponding ChPT ap- 
proximation that, in principle, can be improved order by 
order - eventually including more subtractions. 



II. FITS TO DATA AND LATTICE RESULTS 

As commented before, it has been known for long 
that with the one-loop elastic lAM (the mlAM is almost 
identical) in Eq. (fT6|) it is possible to obtain a remarkable 
description of tttt and Kt: experimental data up to some- 
where below 1 GeV. Simultaneously, the lAM generates 
the poles associated to the /o(600), /9(770), AT* (892) and 
k(800) resonances and this is achieved using parameters 
compatible with those of standard ChPT [12l|. However, 
that description was obtained from a fit to experimen- 
tal data, and therefore it is mostly sensitive to the LECs 
Li, L2, L3 that predominantly govern the s dependence of 
partial waves, but much less so to the rest of LECs that 
carry an explicit meson mass dependence. Of course, 
since now we want to extrapolate the lAM fits to non- 
physical masses, it is very important that we use a good 
description of the mass dependence in observables like 
masses, decay constants, etc. before extracting conclu- 
sions about resonance behavior. For that reason, we are 
presenting here an updated lAM description of exper- 
imental data simultaneously fitted to the available lat- 
tice results on the mass dependence of M^//^, M^//x 
and Mfi / fx as well as scattering lengths for the doubly 
charged channel in tttt, Ktt and KK scattering. Note 
that for the moment, these lattice data are only avail- 



able in the highest isospin combination for each particle 
pair. 

In order to change the masses and decay constants ac- 
cording to Eqs.®,®,® and (HUl), we need first to ex- 
tract the tree level quantities: Mq^, and /o from 
the physical values of the pion, kaon and eta masses as 
well as the three decay constants J-mfx and /,,. Note 
that Mq^ will be obtained from the Gell-Mann-Okubo 



relation: '^M^j^ — Mq^ — 3Mg^ = 0. Since there are more 
physical values than tree level constants, for a given set of 
LECs we actually use the tree level constants that best fit 
the physical ones. Thus, the physical masses and decay 
constants that we will obtain when recovering them from 
the tree level ones will be only approximate. This is, of 
course, the consequence of using a truncated expansion 
- ChPT to one-loop ~ to describe observables. 

We have made two fits whose resulting LECs sets are 
given in the two last columns of Table H] Since there are 
many parameters, there are strong correlations. Thus, 
sets with quite different parameters can give raise to 
acceptable descriptions of data, depending on how one 
weights experiment and lattice results. On Fit I we have 
fitted to experimental data coming from [l^ and to lat- 
tice results given in [l^. We show in Fig.l these data, for 
the [I, J) = (0,0) and (2,0) waves, where many differ- 
ent experiments are actually incompatible, but we have 
fitted to the phase shifts arising from the dispersive anal- 
ysis of the experimental data in [20I] , where a complete 
set of Forward Dispersion Relations and Roy Eqs. was 
constrained on a phenomenological fit to all waves. For 
the (1,1) wave we have used also the phenomenological 
phase shifts from that solution since, apart from the dis- 
persive constraints, it fits the data of the electromagnetic 
form factor of the pion, which is much more reliable and 
precise than the existing experiments on (1,1) pion-pion 
scattering. Anyway, since for ttK and other waves we 
are still using scattering, and also because the method 
has an intrinsic error due to the NLO approximation on 
the integrals, we have added in quadrature to the exper- 
imental data errors a constant error of 2 degrees and a 
variable error of 5% of the phase shift and to lattice re- 
sults on masses over decay constants 5% of their values 
also in quadrature to their errors. We have also intro- 
duced a constraint so that the LECs don't differ much 
from those found in the K14, analysis to two loops of 
by weighting also in the the LECs with the values in 
3. On Fit II we have given an additional weight to the 
large \/Nc constraint 2Li — L2 — (dividing its error by 
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FIG. 2: Results of our lAM fits versus experimental data on tttt and -wK scattering. The continuous and dashed lines 
correspond, respectively, to Fits I and II, whose parameters are given in Tabled For comparison we show the results of the 
lAM if we used the ChPT LECs obtained from the two- loop analysis of Kin, decays listed also in Table U (dotted line) as well 
as the results of standard non-unitarized ChPT with the same set of LECs (dot-dashed line). The plotted data correspond to 
experimental results [T^, which are often incompatible. For that reason, in our fits, and for tttt scattering, we have actually 
used the results of a dispersive analysis of these data [20| . 



10 when calculating the y^) whereas we have relaxed the 
constrains on 5ii and (5i/2o (dividing their by 1.5). 

For comparison, also in Table U we provide three typi- 
cal sets of LECs available in the literature obtained from 
data analyses using dispersive techniques plus ChPT. 
Those on the first and the second columns, come from 
a one and two- loop analysis of Kn decays .8], where 
L4 and Lg were set equal to zero (following leading or- 
der 1/Nc arguments). The "Roy-Steiner" column comes 
from a dispersive analysis of irK scattering [9]. Note 
that the LECs in these sets are frequently within more 
than two standard deviations from one another, and we 
consider that their difference is indicative of the typi- 
cal size of systematic uncertainties in our knowledge of 
LECs. As commented above, since the one-loop lAM 
generates correctly only the s-channel leading logs of the 
two-loop calculation, which are dominant at low ener- 
gies, it is not clear whether we should compare with the 
LECs obtained in the one or two loop ChPT analysis. 
Actually, all of our lAM LECs lie very close, or within 
the uncertainties, of at least one of the previous deter- 
minations given in the Table. Taking into account the 
uncertainties in these non-unitarized determinations, we 
consider that the agreement between the lAM LECs and 
previous determinations is fair. Let us remark that the 



relevant fact about this comparison is to note that we 
do not need to make any fine tuning of the LECs, like 
changing well established signs, changing order of magni- 
tude, etc. , to describe the experimental and lattice data 
simultaneously. 

Finally, we also provide in Table U the lAMIII set of 
LECs, which corresponds to one of the three fits obtained 
using the coupled channel lAM in [l3|. This set was 
fitted to experimental data only and the uncertainties 
quoted are just statistical. Taking into account that we 
are using the single channel lAM instead of the coupled 
channel one, and the estimate of systematic uncertain- 
ties discussed above, we see that our new fits including 
new experimental data and lattice results, are not too 
different from those already obtained in 

hi Fig[2]we show the results of our fits compared with 
experimental data on tttt and i:K elastic scattering phase 
shifts. The best description is given by Fit I (continuous 
line), whereas Fit II gives a somewhat too heavy p(770) 
vector resonance (by roughly 50 MeV, i.e., a 6% error). 
For comparison, we show as a dotted line the results of 
the lAM if we used the ChPT LECs obtained from the 
two-loop analysis of Kn decays listed in Table HI We also 
show as a dot-dashed line the results that would be ob- 
tained if the non-unitarized ChPT one loop results are ex- 



8 



trapolated to higher energies using the same set of LECs. 
Note that the lAM resuhs describe rather weh both the 
resonant and non resonant shapes up to 1 GeV or shghtly 
above, except for the scalar isoscalar Sao, that is only de- 
scribed up to 800 MeV. This is due to the presence of the 
sharp rise caused by the /o(980) resonance, that decays 
mostly to two kaons and can only be described with the 
coupled channel lAM formalism, [10, 12, 16], that we do 
not use here for the reasons explained above. 

Those results are, of course, well known, and these fits 
would just be an update of [lo| if we had not also included 
lattice data on the fit, that we show in FiglH Note that 
we are fitting results on AIt^/ fw, Mt^/ fx and Mk / Jk and 
the 7r+7r+, K^K^ , K^tt^ scattering lengths [l9|. Once 
again we show Fits I and II as continuous and dashed 
lines, respectively, together with lAM results using the 
LECs from the two-loop analysis of Kn decays listed in 
Table U (dotted line) and non-unitarized ChPT to one- 
loop with the same set of LECs (dot-dashed line). As 
explained above, we do not consider that our method 
should be trusted for pion masses heavier than 440 MeV, 
being optimistic, and that is why the heavier mass region 
is shown as a grey area. 

III. DEPENDENCE ON U AND d QUARK 
MASSES 

Now that we have a good description of both the energy 
dependence of pion-pion amplitudes together with the 
mass dependence of the few observables available from 
lattice, we can change the value of the light quark mass, 
keeping fixed, and predict the behavior of the reso- 
nances generated within the lAM. 

A. Light vector mesons: the p(770) and /S'*(892) 

The /9(770) and AT* (892) vector resonances are well es- 
tablished qq states belonging to an SU(3) octet. The first 
is produced in tttt scattering, and its quark mass depen- 
dence was aheady studied within SU(2) ChPT 0. Here 
we will just check that we reobtain very similar results 
within the SU(3) formalism, while describing simultane- 
ously the lattice observables shown in FiglH However, 
the A'* (892) appears in ttK scattering and can only be 
obtained using SU(3) ChPT as we do here. 

1. Mass and width 

Thus, in FigHJ we show the dependence of the light 
vector resonances on the non-strange quark masses, us- 
ing one-loop SU(3) ChPT unitarized with the lAM. For 
each resonance, these masses and widths are defined from 
the position of their associated pole in the second Rie- 
mann sheet, through the usual Breit-Wigner identifica- 
tion: yfSpoie = M — i r/2. We show the results for Fits I 



and II as continuous and dashed lines, respectively. The 
results for both fits are very consistent and their differ- 
ence can be taken as an estimation for systematic uncer- 
tainties in the choice of LECs. To suppress systematic 
uncertainties we give all quantities normalized to their 
physical values. Note that we provide two scales for the 
mass variation: In the upper horizontal axis, we show 
the variation of the quark mass in terms of ih/wP^^^, 
whereas in the lower horizontal axis we show the varia- 
tion of the pion mass in terms of Mj^/Mp^^^. The one- 
loop ChPT relation between these two scales is given by 
Eqs.([2]) and ([3]). To be precise, this relation changes for 
different LECs, but, as we already showed in FiglTJ the 
difference is too small to be observed with the naked eye 
in the axes of Fig|31 

In the left panels we also show, as a dotted line, the 
SU(2) ChPT result aheady obtained in [|, which, is 
fairly consistent with the new SU(3) results. Of course, 
the difference is somewhat larger when the pion mass 
is closer to the kaon mass, and the kaons start playing 
a more prominent role. Of course, since the SU(2) re- 
sults [3,] already described fairly well the available lattice 
calculations for the /o(770) mass, so it happens with the 
SU(3) results here. In addition, this ensures that the Mp 
dependence on Mjr agrees nicely with the estimations for 
the two first coefficients of its chiral expansion [2l[ , which 
was already checked in the SU(2) case Q. 

Since the vertical scale is the same for the p{77Q) and 
A'* (892) plots, the similarity of their behavior is very ev- 
ident. Both their masses increase smoothly as the quark 
mass increases, but much slower than the pion mass. 
Some differences can be observed for small m, but this 
is due to the fact that the SU(3) breaking between the 
p{770) and the AT* (892) is more evident since we keep rUg 
fixed to its large physical value. What is interesting to 
observe is that the naive rule of thumb frequently used in 
the literature that dMn/dm = A^]^, where N]^ is the 
number of valence non-strange quarks, yields the correct 
order of magnitude, (and this is how it has been used in 
22]) but would predict a 2:1 relation for the slope of the 
p(770) with respect to that of the A'* (892), which is not 
observed for light quarks. 

Continuing with our analysis, we note that, as the 
quark mass increases, the two-pion and pion-kaon thresh- 
old grow faster than the masses of the resonances and, 
as a consequence, there is a strong phase space suppres- 
sion than can account exclusively for the decrease of their 
widths. We show in the lower panels the dependence 
of T p and P^. normalized to their physical values. The 
decrease of the widths is largely kinematical, following re- 
markably well the expected reduction from phase space 
as the masses of the NGB increase (thin continuous and 
dashed lines corresponding to Fit I and Fit II respec- 
tively, although for the p(770) they overlap so well that 
the thin lines are not seen ). This result was already 
found for the p{77Q) within the SU(2) formalism and is 
nicely confirmed here. This suggests that there is no dy- 
namical effect through the vector coupling to two mesons. 




FIG. 3: Result of the unitarized fits to lattice calculations of Mt^I J-k^ Mt^/ Jk, Mk / Jk and the ■k'^-k^ , K'^K^ , K'^n'^ scattering 
lengths. The continuous and dashed lines correspond, respectively, to Fits I and II, whose parameters are given in Tabled For 
comparison we show the results of the lAM if we used the ChPT LECs obtained from the two-loop analysis of Ku decays listed 
also in Table |l] (dotted line) as well as the results of standard non-unitarized ChPT with the same set of LECs (dot-dashed 
line). Lattice results come from [l^. The grey area lies beyond our applicability region, however, it is useful to check that our 
description does not deteriorate too rapidly. 



as we will analyze next. 



2. Coupling to two mesons 

The dynamics of resonance-meson-meson interaction is 
encoded in the coupling constant that we obtain from the 
residue of the amplitude at the pole position as follows; 
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5polc, 



t(s) 



(19) 



where the normalization factors are chosen to recover the 
usual expression for the two-meson width of narrow vec- 
tor resonances: 



|5P 



(20) 



|k| being the modulus of the meson three- momentum. 
Actually, by identifying ^Spoio = My — iFy, we have 
explicitly checked that we obtain the same numerical 
value for the coupling with both equations. We find 
IOp-rttI - 6.1 and \gK'-^K\ - 5.5. 

Then, on FiglS] we show the dependence of the ^ptttt 
(left) and the qk'ttK (right) couplings with respect to 



the pion mass (lower horizontal scale) or the non-strange 
quark mass m (upper horizontal scale) . In order to sup- 
press systematic uncertainties, we have normalized the 
couplings to their physical values. Note that the gp^ir 
is remarkably constant, deviating from its physical value 
by 2% at most, despite the fact that the quark mass is 
changed by a factor of 9. It is also relevant because it jus- 
tifies the constancy a.ssumption made in lattice studies of 
the p(770) width [23|. The gK'-r^K is also quite indepen- 
dent of the non-strange quark mass, deviating by 10% at 
most in the chiral limit and by less than 4% when the 
quark mass is increased by a factor of nine. The results 
for Fit I and II are almost indistinguishable. 

The constancy of the vector-meson-meson couplings, 
together with the classic KSRF relation [24,], provides a 
striking connection between the quark mass dependence 
of the rho mass and the pion decay constant. Actually, 
the KSRF relation, obtained from PCAC and vector me- 
son dominance, reads: 



(21) 



Note that in our calculation we are obtaining Mp from a 
one-loop ChPT unitarized calculation, whereas comes 
simply from the next to leading order ChPT calculation. 



10 



A A 

m / m 



phys 

0.5 1 2 3456789 




1 1.5 2 



A A 

m / m 



0.5 1 



phys 
2 3 4 5 



6 7 8 9 



1 1 1 r 




0.5 



1 1.5 2 



2.5 



1.4 



0.5 1 



A A 

m / m 



phys 

2 3 4 5 6 7 



8 9 



T r 



"1 1 1 1 — I — r 




1 1.5 2 



1.4 



A A 

m / m 



0.5 1 



phys 
2 3 4 5 



6 7 8 9 



"1 1 1 1 1 r 



T — r 




1 1.5 2 



FIG. 4: Dependence of the p(770) and if* (892) mass and width with respect to the non-strange quark mass m (horizontal 
upper scale), or the pion mass (horizontal lower scale). Note that we give all quantities normalized to their physical values. 
The thick continuous and dashed lines correspond, respectively, to Fit I and Fit II described in the text with unitarized SU(3) 
ChPT. For the p these results are very compatible with those in [3] using SU(2) ChPT (dotted line). The continuous (dashed) 
thin line shows the Mtt dependence of the widths from the change of phase space only, assuming a constant coupling of the 
resonances to two mesons, p(770) to tttt and A"* (892) to irK, calculated from the dependence of masses and momenta given by 
Fit I (II). For the p(770) the thin and dashed lines overlap completely. 



but, of course, without unitarization. It is therefore quite 
remarkable that the ratio Mp/ f^^ obtained from our am- 
plitudes, shown in FiglHl is constant within less than 5% 
accuracy, when the quark mass varies by a factor of 9, or 
the pion mass by a factor of 3. Note that, as usual, in 
Figiniwe have normalized the ratio to its physical value. 
It seems that the simple KSRF relation holds remarkably 
well up to surprisingly large values of the non-strange 
quark mass, and therefore the Alp quark mass depen- 
dence can be recast with the same factor as that for 

A similar result is found for the if* (892) whose ra- 
tio Mk* //tt is also shown in FiglB] to deviate by less 



than 2% from its physical value. Note that, accord- 
ing to the second reference in '24'], the fx dependence 
does not show up in the relation. Actually, had we used 
Mk" I V/Tr/if instead, the deviation would have been a 
factor of 3 larger. 



B. Light scalar mesons: the /o(600) and k(800) 

The /o(600), or sigma, and the k(800) scalar mesons 
are still somewhat controversial. The main problem is 
their huge width that makes their experimental identi- 
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FIG. 5: Two-meson- vector coupling dependence with respect to the non-strange quark mass m (horizontal upper scale), or the 
pion mass (horizontal lower scale). Note we normalize the couplings to their physical values. We show on the left the p(770) 
coupling to two pions and on the right that of the K* (892) to ttK (continuous and dashed lines correspond to Fit I and Fit II 
respectively). 
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FIG. 6: Ratio of vector resonance masses to the pion decay constant dependence on the non-strange quark mass (horizontal 
upper scale), or the pion mass (horizontal lower scale). Once again, we normalize all quantities to their physical values. We 
show the Mp//^ on the left and Mk* / f-n on the right (continuous and dashed lines correspond to Fit I and Fit II respectively). 
Both masses seem to follow the quark mass dependence up to less than 5% and 2%, respectively. 



fication complicated. Despite the fact that their pole 
mass and width has been determined by several groups 
vifith the help of model independent dispersive techniques 
(with and without ChPT input) and a fairly reasonable 
agreement (see [l,[i3,[H| for recent determinations), they 
are still cited with extremely cautious and conservative 
estimates in the PDG 26]. Their nature is even more 
controversial, and as commented above, there are no 
present lattice calculations with realistic quark masses 
that could shed some light on the problem. It is therefore 



even more interesting to obtain predictions on their quark 
mass dependence. Compared with the vector case, there 
is an additional complication because now we do not nec- 
essarily expect a similar behavior between the «;(800) and 
the /o(600), since although the former should belong to 
an SU(3) octet, the latter could be in the singlet, the 
octet, or have a significant mixture of both. As a matter 
of fact, there are indications that its singlet component 
is actually dominant [13, HI] . 
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1. Mass and width 

As the data shows in Fig. [2l the sigma and kappa 
resonances do not present a peak nor a Breit-Wigner 
shape in the meson-meson scattering (/, J) = (0, 0) and 
(1/2,0) waves, respectively. Once again, these masses 
and widths are defined from the position of their as- 
sociated pole in the second Riemann sheet, as follows: 
y/Spoie = Ajf — ir/2, but one should keep in mind 
that these scalar states do not present the typical Breit 
Wigner shape, so there is no immediate equivalence of 
the mass in terms of a peak in the cross section or a time 
delay in the propagation. 

In Figl7]we show the pole mass and width dependence 
of light scalar resonances on the non-strange quark mass. 
As in FiglH we show quantities normalized to their phys- 
ical values and we provide two scales for the horizontal 
axis: m/mphys (upper horizontal axis) and MT^/MT^phys 
(lower horizontal axis). Once again, the continuous line 
represents the results for Fit I, the dashed line those of 
Fit II, and the dotted line stands for the results of unita- 
rized SU(2) ChPT for the /o(600). As before we find that 
the Fit I and II are very consistent with each other, and, 
for the /o(600) also with the existing SU(2) calculation 
of i. 

The most prominent feature of the scalars behavior 
is the appearance of two branches for the mass as de- 
fined above, already observed for the cr in 0]. The rea- 
son is that for physical values of the quark mass, the 
poles associated with resonances appear as conjugated 
poles in the second Riemann sheet, i.e., there are poles 
at y/Spoie = MzLi r/2. Of course, only the one in the low- 
est half plane is continuous with the physical amplitude 
in the real axis, and this is the one responsible for the 
physical resonance. However, as the quark mass increases 
these poles move closer to the real axis until they join in 
a single pole below threshold, but still in the second Rie- 
mann sheet. If the quark mass is increased further, the 
poles split again but without leaving the real axis. The 
position of each one of these poles corresponds to each 
one of the branches that we show in the upper panel of 
FigEl 

Although this qualitative behavior is a well known pos- 
sibility for potentials in scalar channels, one-loop unita- 
rized ChPT is predicting the quark mass value for which 
it occurs, which is a genuine prediction for QCD. For 
scalar-isoscalar tttt scattering it was already observed in 
the SU(2) case @. Here we are confirming this posi- 
tion when using SU(3) instead of SU(2) ChPT, but we 
see it also happening for the k(800), although the point 
at which it happens depends more on the set of LECs. 
For this reason, we think that the existence of this non- 
analyticity of the k(800) pole is robust, but not so much 
the precise quark mass value where it occurs. 

This "apparent splitting" cannot occur for higher par- 
tial waves since they all carry a k^"^ factor that forces the 
conjugated poles to join the real axis exactly at thresh- 
old, and then one of them jumps to the first Riemann 



sheet. 

Apart from the evident qualitative similarities between 
the behavior of the /o (600) and the kappa, it is also clear 
that quantitatively they behave somewhat differently. In 
particular, the growth of the k(800) mass before the 
"splitting point" is much softer than for the /o(600), and 
even softer than the p{770) and if* (892) growth shown 
in FigU (please note the difference in scales between both 
figures) . 

In the lower panels of Fig|7] we show the quark mass 
dependence of the sigma and kappa widths. On the left 
we show that the decrease of the sigma width we find with 
the SU(3) one-loop lAM is very consistent between Fits I 
and II, and confirm the previous results within SU(2) Q. 
On the right we show the results for the k(800) width. 
We also show that the width decrease for both of them 
cannot be attributed to the phase space reduction, due 
to the increase of pion and kaon masses, naively expected 
from the narrow width approximation 

which we show as a thin continuous (dashed) line cor- 
responding to Fit I (II). Although the shape of the de- 
crease is slightly different for the a and k, both scalars 
behave very differently than vector mesons. Actually, we 
will see next that this implies that the scalar couplings to 
two mesons have a much stronger quark mass dependence 
than the vector ones. 



8. Coupling to two mesons 

As we have just seen, the narrow width approximation 
in Eq. ([22l) above is of little use for scalars. But, of course, 
we can still extract the coupling constant from the residue 
as we did for vectors, although now the equation reads 

= -167r lim (s - Spoic) t{s), (23) 

S^Spolo 

We find \g^^^\ ~ 2.86 GeV and \g^^K\ ^ 3.6 GeV, to 
be compared to Igcnr-jrl — 2.97 ±0.04 GeV and I^KTr/fl — 
4.94±0.07 GeV, obtained in [Hi or the ~ 2.2 aver- 

age obtained in [2^. The agreement is fairly reasonable, 
taking into account that the data that have been used, 
the a and n poles, and the models in those references 
differ substantially for each reference. 

Thus, in Fig|8] we show the quark mass dependence 
(upper horizontal scale) or pion mass dependence (lower 
horizontal scale) of gaTnr and gmrK- As usual, all quan- 
tities are normalized to their physical values. Compared 
with FiglS] (note the different scales), we see that these 
couplings show a much stronger quark mass dependence. 
Moreover, they increase dramatically near the point of 
the "apparent splitting". Beyond that point there are 
two non-conjugate poles lying on the real axis below 
threshold in the second Riemann sheet. For this rea- 
son, after the splitting point we plot two curves for each 
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FIG. 7: Dependence of the /o(600) and k(800) mass and width with respect to the non-strange quark mass m (horizontal 
upper scale), or the pion mass (horizontal lower scale). Note that we give all quantities normalized to their physical values. 
The thick continuous and dashed lines correspond, respectively, to Fit I and Fit II described in the text with unitarized SU(3) 
ChPT. For the /o(600) these results are very compatible with those in 0] using SU(2) ChPT (dotted line). Let us remark that 
both resonances, being scalar, develop two poles on the real axis for sufficiently high masses. The thin lines show the decrease 
of the widths if it were only due to phase space reduction (the thin continuous line corresponds to Fit I and the dashed one to 
Fit II). 



fit. The lowest curve corresponds to the pole closest to 
the threshold, that eventually jumps into the first Rie- 
mann sheet. This threshold crossing from one sheet to 
the other corresponds to the point where the coupling 
tends to zero in the figures, in good agreement with the 
well known result in Actually this can be checked 

numerically, because, as shown in [31[ the coupling is in- 
versely proportional to the energy derivative of the one- 
loop function ( G(s) in ^ and J(s) in ChPT [12]), which 
is divergent at threshold. Despite this consistency check, 
within our approach this occurs at pion masses close to 
the naive applicability limit, and therefore the exact M^r 



value when this happens is not very reliable. 



IV. DEPENDENCE ON THE STRANGE 
QUARK MASS 

Up to here we have only been changing the values of 
the non-strange quark mass keeping mj fixed. However, 
since we are dealing with the full SU(3) ChPT formalism, 
we are now able to change the strange quark, keeping m 
fixed. The dependence of hadronic observables on the 
strange quark mass is also of interest for lattice studies 
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FIG. 8: Two-meson-scalar coupling dependence with respect to the non-strange quark mass m (horizontal upper scale), or the 
pion mass (horizontal lower scale). Note we normalize the couplings to their physical values. We show on the left the /o(600) 
or "sigma" coupling to two pions and on the right that of the k(800) to irK. Let us remark that the two poles in the real axis 
show different couplings, which explains the doubling of the lines above the double branch point. 



and for cosmological considerations [6|. As we explained 
in Sect ll Al we will only vary the strange quark mass in 
the limited range 0.7 < ms/nis phys < 1-3 to ensm'c that 
the kaon does not become too heavy to spoil the ChPT 
convergence nor too hght to require a coupled channel 
formalism to deal with the K*{892) or k(800) resonances, 
thus introducing additional model dependences in our 
approach. 



A. Light vector mesons: the p(770) and K*{892) 
1. Mass and width 

As in previous sections we define the mass and width 
of the vector resonances from the position of their as- 
sociated poles. Thus, in the upper panels of Figl9] we 
show the quark mass dependence (or kaon mass depen- 
dence in the lower horizontal scale) of the p and i4r*(892) 
masses. In the lower panels we show the dependence of 
their widths. As usual, all quantities are normalized to 
their physical values to suppress systematic uncertainties. 

As it could be expected, both the mass and width of the 
/9(770), being non-strange, are almost independent of the 
strange quark mass within the range of study. Note that 
the p mass actually decreases very slightly, by roughly 
1%. Since the pion mass almost remains constant - see 
Eq.(I31) and the Lg, -^4 values in Table|T]~, this implies that 
phase space decreases slightly for smaller strange quark 
mass and the /5(770) width decreases accordingly. Actu- 
ally, we can check in FiglS]that the width reduction fol- 
lows remarkably well the phase space reduction expected 
from Eg. d^ni) (thin continuous and dashed lines). 



Looking now at the right panels of FigEl we notice 
that, as expected, the X*(892) shows a much stronger 
dependence than the p(770) on the strange quark or the 
kaon masses. On the one hand, when the kaon mass 
is made lighter, the ii'*(892) mass decreases, as it hap- 
pened when changing the light quark mass, although 
much faster, i.e., up to 5% when the kaon mass decreases 
by 20%. Nevertheless, and contrary to what happened 
when reducing m, the i^T* (892) width increases signif- 
icantly, up to 40%. This is due to the fact that the 
K*{8%2) decays to ttX, but the kaon mass decrease is 
faster than that of the if* (892). On the other hand, when 
the kaon mass is made heavier, the if* (892) mass grows, 
but much slower than the kaon mass, so that phase space 
shrinks and the resonance width decreases once more. 
We are also showing as thin lines the expected variation 
of the widths if their only quark mass dependence came 
from the change in the particles masses and the naive 
phase space suppression in Eg. ([20]) (thin continuous for 
Fit I and thin dashed for Fit II). We see that they are 
in very good agreement with our results from the lAM, 
which suggests that their coupling to two mesons is al- 
most independent of the quark masses, which we will see 
next. 



2. Coupling to two mesons 

Thus, in Fig. [TU] we show the dependence both on rria 
and kaon masses of the vector to meson-meson couplings. 
As usual everything is normalized to their physical val- 
ues. It can be noted that within the range of variation 
under study, which is 30% for the strange quark mass 
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FIG. 9: Dependence of the p(770) and K*{W2) mass and width with respect to the strange quark mass rus (horizontal upper 
scale), or the kaon mass (horizontal lower scale). Note that we give all quantities normalized to their physical values. The thick 
continuous and dashed lines correspond, respectively, to Fit I and Fit II described in the text with unitarized SU(3) ChPT. The 
thin lines show the decrease of the widths as if it were only due to phase space reduction (the thin continuous line corresponds 
to Fit I and the dashed one to Fit II). 



in either direction, both the gp^-n and gK'-nK couplings 
change by 1% at most. 

In Fig. [TT] we show the results for the KSRF relation 
variation in terms of the strange quark mass. Since the p 
coupling has virtually no dependence on m^, the relation 
remains trivially constant. For the i4r*(892) the relation 
is well satisfied (to within less than 5% from the physical 
value) in the whole rn^ range of our study. 

B. Light scalar mesons: the /o(600) and k(800) 

We simply repeat the procedure we used to study the 
light quark variation in Sec. IIII Bl but this time changing 



the strange quark mass instead, and keeping rh fixed. 



1. Mass and width 

Thus, in Fig. 1121 we show the variation of the sigma 
and k(800) masses and widths with respect to the kaon 
mass variation (lower horizontal scale) or the strange 
quark mass (upper horizontal scale). Once again all 
masses are normalized to their physical values. As it 
could be expected, we see in the left panels that the 
change on the sigma is smaller than 1% on both mass 
and width (beware we have changed the scale with re- 
spect to the previous Fig. [7] to make the changes more 
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FIG. 10: Two- meson- vector coupling dependence with respect to the strange quark mass nis (horizontal upper scale), or the 
kaon mass (horizontal lower scale). Note we normalize the couplings to their physical values. We show on the left the p(770) 
coupling to two pions and on the right that of the K* (892) to ivK. 
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FIG. 11: Ratio of vector resonance masses to the pion decay constant dependence on the strange quark mass (horizontal upper 
scale), or the kaon mass (horizontal lower scale). Once again, we normalize all quantities to their physical values. We show the 

Mp/fTT on the left and Mk* / fir on the right (continuous and dashed lines correspond to Fit I and Fit II respectively). For 
the p the result is trivial since we have already shown that its coupling does not depend on rUa- However the _ftr*(892) mass 
deviates from the /,r quark mass dependence by less than 5% with respect to the physical value. 



visible). 

A much bigger effect is seen for the /t(800) in the right 
panels, whose mass changes by as much as 12% from 
its physical value within the range of study, whereas the 
width changes by as much as 20%. However, its mass 
dependence, despite being somewhat stronger than for 
its vector counterpart K*{9,92), is still softer than for the 
kaon itself. This is the reason why, as the k;(800) becomes 
lighter its width increases, and viceversa. 



In the lower panels we have also plotted the expected 
naive phase space reduction. This time, however, as 
the sigma properties barely depend on the strange quark 
mass, we only sec a significant deviation from that naive 
behavior in the case of the k(800). 
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FIG. 12: Dependence of the /o(600) and ft:(800) mass and width with respect to the strange quark mass (horizontal upper 
scale), or the kaon mass (horizontal lower scale). Note that we give all quantities normalized to their physical values. The thick 
continuous and dashed lines correspond, respectively, to Fit I and Fit II described in the text with unitarized SU(3) ChPT. The 
thin lines show the decrease of the widths as if it were only due to phase space reduction (the thin continuous line corresponds 
to Fit I and the dashed one to Fit II) 



2. Coupling to two mesons 



V. SUMMARY AND CONCLUSIONS 



For all means and purposes, with respect to strange 
quark mass variations, the sigma coupling to two mesons 
turns out to be a constant within our approximation, as 
can be seen in the left panel of Fig. [131 



In contrast, the g^-nK coupling shows some dependence 
on the strange quark mass. Actually, it grows by 6% 
when the kaon mass is increased by 18% from its physical 
value. 



In this work we have studied the quark mass depen- 
dence of the light vector and scalar resonances generated 
as poles of meson-meson scattering elastic amplitudes 
within unitarized one-loop Chiral Perturbation Theory 
(ChPT). This dependence is of interest to relate lattice 
results to hadronic observables, but also for anthropic 
and cosmological considerations. The use of an SU(3) 
formalism extends previous studies within SU(2), allow- 
ing us to study the behavior of strange resonances like 
the k(800) and if* (892), but also to study variations not 
only of the light u and d quark masses, but also of the 
strange quark mass. 
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FIG. 13: Two-meson-scalar coupling dependence with respect to the strange quark mass rua (horizontal upper scale), or the 
kaon mass (horizontal lower scale). Note we normalize the couplings to their physical values. We show on the left the /o(600) 
or "sigma" coupling to two pions and on the right that of the k;(800) to nK. 



After a brief introduction on how ChPT provides a 
model independent expansion of pion, kaon and eta 
masses and decay constants, as well as their two body 
interaction amplitudes, we have reviewed how this series 
can be used inside a dispersion theory formalism to con- 
struct the so called Inverse Amplitude Method (lAM) 
amplitudes that satisfy elastic unitarity while respecting 
the ChPT expansion. It has been known for long that 
the elastic lAM reproduces well the meson-meson elas- 
tic scattering data up to 800-1000 MeV, including the 
resonance region. Note that we have refrained for the 
moment to use the very successful coupled channel lAM 
precisely because at present it lacks a dispersive deriva- 
tion, and we want to avoid as much model dependence 
as possible. Of course the experimental data may fix 
rather well the energy dependence but not so well the 
mass dependence. For that reason we have presented here 
a new lAM analysis including simultaneously the exist- 
ing lattice results on meson masses, decay constants and 
scattering lengths. We obtain a fairly good description 
of experiment and lattice data using chiral parameters 
rather similar to existing one and two-loop determina- 
tions. No fine tuning of parameters is required. Once 
this is done, we have varied the quark masses within cer- 
tain ranges that ensure the applicability of the elastic 
lAM for the resonances under study: m/wphys < 9 and 
0.7 < ms/nis phys < 1.3 (m is the average mass of the u 
and d quarks). In practice, in ChPT we have changed 
the squared pion and kaon masses, which, at leading 
order, are proportional to quark masses. Although we 
have shown in Fig[T]that this simple approximation works 
within roughly 10% accuracy, we have carefully included 
the full one-loop corrections, and shown the quark and 
meson mass variation independently in all plots. 



In the second Riemann sheet of these amplitudes, the 
lAM generates the - conjugated pairs of - poles asso- 
ciated to the vector p(770), K*{892) and scalar /o(600) 
and k(800) resonances. Light vector resonances are well 
established and there is little relevance on whether we re- 
fer to their "pole" or Breit-Wigner mass and widths. In 
contrast, the scalar /o(600), or "sigma" and the k(800) 
are rather controversial due to their large apparent width 
and the lack of a Breit-Wigner shape in the meson-meson 
scattering phase shifts. To avoid complications, we have 
always presented our results in terms of "pole" definitions 
of masses, widths and couplings. 

For the /o(600) and p(770) resonances, which appear in 
TTTT scattering, we have nicely confirmed the similar uni- 
tarized one- loop SU(2) ChPT analysis performed in 0]. 
When increasing m both the sigma and p masses grow 
faster than the pion mass, whereas their widths decrease. 
However, the p(770) mass behaves smoothly in the whole 
quark mass range, whereas, roughly at Mtt ~ 340 MeV, 
the /o(600) pole and its conjugated pair meet in the sec- 
ond Riemann sheet below threshold, producing a non 
analyticity - or "apparent splitting" in two branches - 
of the sigma mass in terms of M^. In addition, we con- 
firm that the p(770) width decrease, as rh grows, fol- 
lows remarkably well the simple expectations of phase 
space reduction already found within the SU(2) formal- 
ism. Once again, such a simple behavior is not observed 
for the sigma. 

Of course, the SU(3) formalism allows us now to study 
also the K*{892) and «:(800) resonances in nK scattering. 
We find that both the mass and width of the K*{892) 
behave qualitatively and quantitatively in a very similar 
way to those of the p(770), which could be expected given 
the fact that they belong to the same octet. In addition, 
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we have explicitly calculated here their couplings to two 
mesons, from the residue of the partial wave at their as- 
sociated pole, finding that they are both remarkably in- 
dependent of the non-strange quark mass, as suggested 
from the width behavior. The K*{892) coupling is quite 
well approximated by a constant, although not so well as 
in the p case. This could be of relevance when computing 
its width on the lattice as it has already been done for 
the p m. 

It therefore seems that light quark masses play no 
significant role in the dynamics of the dominant decay 
modes of vector mesons, namely p — >■ tttt and K* — > irK, 
since their couplings seem to be independent of light 
quark masses and all their width variation can be at- 
tributed to the phase space modification due to changes 
in the masses of all particles. 

Furthermore, this provides a hint, checked here by ex- 
plicit calculation, that the KSRF relation, that approx- 
imates these couplings by 5 ~ Afv/2\/2/7r, holds to less 
than 5% when changing m from to 9 times its physi- 
cal value. It is remarkable that this relation is so well 
satisfied, first, because ours is a one-loop calculation, 
which, in principle includes higher order pion mass cor- 
rections to KSRF, and the pion mass becomes rather 
large, but, second, because our resonance masses come 
from unitarized amplitudes whereas /^r stems from the 
non-unitarized ChPT truncated series. 

Concerning the k(800), its behavior is qualitatively 
similar to that of the "sigma" , including the "apparent 
mass splitting" in two branches, which is a feature that 
can only occur for scalars. However, the k(800) non- 
strange quark mass dependence is softer than for the 
sigma. Still the pion mass where the k(800) "apparent 
mass splitting" occurs is similar to that of the sigma, al- 
though with bigger uncertainties ^ 340 — 400 MeV. 
Of course, contrary to the vector case, one could now ex- 
pect some differences between the two scalars since they 
do not necessarily belong to the same octet and actu- 
ally, the sigma is believed to be predominantly the singlet 
state [13, , and it could even allow for a glueball com- 
ponent. As we did with the vectors, in this work we have 
also calculated explicitly the behavior of the scalar cou- 
plings to two mesons under quark mass variations. We 
find a qualitatively similar behavior for both gan-K and 
gKTvK- contrary to vectors, they cannot be considered 
constant within the variation range, particularly when 
comes close to the "apparent mass splitting" value, 
where it suffers a dramatic enhancement. 

Finally, since we use the SU(3) formalism, we have 
been able to study the dependence of light resonance 
properties on the strange quark mass. Due to the fact 
that the physical mass of the kaons is already quite high 



but also because we want the Mk + M^i threshold to be 
significantly above the K*{892) mass, we have limited 
our study to the range 0.7 < ms/nisphys < 1-3. As it 
could be naively expected, and in contrast to strange res- 
onances, the masses and widths of both the non-strange 
p and a are remarkably independent of the strange quark 
mass. This time, the k(800) mass has a much stronger 
dependence than that of the K*{892) - actually, it grows 
a factor of three faster. Once again, the i^*(892) width 
follows remarkably well the behavior dictated by phase 
space only, and we have checked that its ttK coupling is 
almost independent of mg. The KSRF relation is also 
a fairly good approximation in the whole energy range, 
although not as good as in the case of the non-strange 
quark. Concerning the k(800), once again its coupling is 
strongly dependent on the quark mass, so that its width 
does not follow the naive phase space behavior. 

In summary, we have presented an exhaustive study 
on the strange and non-strange quark mass dependence 
of light scalar and vector resonances appearing in elastic 
Goldstone bosons scattering. For the future, this work 
could be extended to other light scalar mesons like the 
/o(980) and ao(980) using a coupled channel formalism, 
that is somewhat less rigorous as it has no dispersive 
derivation, and also is much more complicated due to 
the presence of the KK threshold. 

To conclude, and apart from the interest for studies of 
constraints on hadronic properties from cosmological or 
anthropic considerations, we think that the quark mass 
dependence studied here will be within the reach of lat- 
tice studies in the not too distant future - it is already so 
for the p meson ~ and we expect our results to be useful 
in the chiral extrapolation of lattice results to physical 
values. 
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